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1 Introduction 



Holography relates strongly-coupled gauge theories to weakly-coupled theories in higher 
dimensions. The original conjecture [1] connects type IIB string theory in the gravita- 
tional AdSs x S5 background to the 4D Af = 4 SU (N c ) superconformal field theory This 
correspondence can be extended to other asymptotically-AdS [2] spaces, and examples 
of geometries in which conformal invariance is broken in the IR are known (see, for in- 
stance, [3] and references therein). From a phenomenological point of view, it is often 
sufficient to consider the simpler relation between non-supersymmetric 5D field theories in 
a Randall-Sundrum [4] background and 4D large-iV c strongly-coupled theories with con- 
formal invariance spontaneously broken [5]. This can be applied, for example, in studying 
properties of QCD at large N c [6] or to clarify the physics of electroweak symmetry breaking 
by strong dynamics [7]. 

The 5D setup involves a slice of AdSs truncated by two branes at z — z\jy and z = 
Zir. The bulk hosts 5D gauge and, possibly, other spin fields, which represent composite 
operators of the dual CFT. One way to probe the dynamics of such theories is to compute 
UV boundary correlators of the 5D bulk fields. For example, the 1PI two-point correlation 
function of boundary gauge fields is given by the expression 



J d*xe-^(A»(x)A,(o)) 1Pl = f- Via> + P -^A u( P 2 ), n(p) - dzG{z ' p) 



G(z,p) 



, (1.1) 

=zuv 



where G(z,p) is a solution of the 5D equations of motion subject to appropriate boundary 
conditions at the IR brane. The holographic dictionary relates this boundary correlator 
to the connected two-point correlation function of a conserved global symmetry current of 
the 4D CFT. The poles of this correlation function are interpreted as resonances of the 
CFT. 

There exists another, seemingly unrelated approach to computing correlation functions 
in 4D strongly coupled large- N c theories, which was proposed long ago by Migdal [8]. 
The program of Migdal aims at reproducing the gauge theory correlators at low p 2 using 
information about the deep Euclidean regime. The main input is the non-analytic be- 
haviour of the correlation functions at large Euclidean momenta, where the correlators, 
at leading order, exhibit a conformal behaviour, lim p 2_ > _ 00 Tl(p) = f{p 2 ) ~ — p 2n log(— p 2 ). 
This asymptotic expression is approximated by a ratio of two polynomials of degree N, 
f(p 2 ) ~ Rn(p 2 )/ Sn(p 2 ), by means of a Pade approximation. Finally, the result is analyti- 
cally continued to low time-like p 2 and the large- N limit is taken. 

It was recently pointed out by Shifman [9] and Voloshin, and shown in detail by Er- 
lich et al. [10], that Migdal's procedure gives results similar to those of 5D holographic 
computations. Indeed, performing a series of Pade approximations on the input function 
f{p 2 ) = — log(— p 2 /fi 2 ), one obtains [10], for large N, 

f( 2, -\og(p 2 /^ 2 )J (2Np/ f ,)+nY (2Np/ f ,) 

The same result would be obtained if we performed a computation of the boundary gauge 
field correlator in the Randall-Sundrum spacetime, with z IR = 2N/fi and z\jy = l//x — > 0. 
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In fact, the similarity of the two approaches is not restricted to A0IS5. For example, 
Pade approximations of the input function f(p 2 ) = (— p 2 //i 2 ) -1 / 2 lead to a result that 
again coincides with the 5D boundary gauge correlator, but this time computed in the 5D 
Minkowski background. 1 

This coincidence seems quite mysterious, as the Migdal approach never makes any 
reference to extra dimensions. Of course, the main ingredient of both methods, and the 
one that allows us to make contact with a large- N c theory, is that the correlation functions 
are meromorphic. But this does not determine them uniquely. Actually, the success of the 
Migdal approach and its relation to holographic calculations raises a number of questions: 

1. Why, in the first place, can the series of Pade approximants be interpreted as physical 
correlation functions of a large- N c gauge theory? The physical interpretation of the 
results obtained by Migdal is possible because the Pade approximants have only 
simple poles, with negative residues on the positive real p 2 axis. This is not a generic 
feature of Pade approximants. 

2. Why, in the large- N limit, do we recover correlators characteristic of theories in 5D? 
The distinguishing feature of 5D models is locality in the fifth coordinate. How is 
this locality encoded in the Migdal approach? 

3. Which 5D geometries and which IR boundary conditions can be reproduced by the 
Migdal approach? 

4. Does the Migdal procedure at finite N also correspond to some physical setup? There 
is a tantalizing similarity between the correlators computed in the Migdal approach 
and those obtained in deconstruction [12]. Both express the polarization function 
as a ratio of two polynomials in p 2 , which converges to a non-analytic function for 
large N and large Euclidean momentum (see [13] for a computation in deconstructed 
AdSs). Is there a precise quantitative connection between the two approaches? 

In this paper we provide answers to these questions. The first one was, in fact, already 
addressed in Migdal's original papers [8]. It turns out that the nice physical properties of 
Pade approximants arise because the input function belongs to the class of so-called Stieltjes 
functions. In the mathematical literature, Pade approximants of Stieltjes functions have 
been extensively studied, and a connection to orthogonal polynomials has been established. 
The physical properties of Pade approximants are intimately related to certain familiar 
properties of orthogonal polynomials. 

One well-known property of orthogonal polynomials is that they satisfy a second-order 
recursion relation in the polynomial degree: 7Tjv+i(w) = (a n w + b n )TT^{w) — c n rc^^i{w). 
We show that the connection between Pade approximation and orthogonal polynomials 
implies that the polynomials entering the Pade approximants also satisfy a second-order 
recursion relation: 

T N+l {p 2 ) = A N {p 2 )T N {p 2 ) - B N {p 2 )T N ^{p 2 ) , (1.3) 
1 This was first observed by Gherghetta, Pomarol and Rattazzi [11]. 
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where T = R,S. Hence, Pade approximants implement automatically a form of locality in 
the discrete space of the polynomial degrees. For large N, this recursion relation reduces to 
a second order differential equation, whose form is analogous to equations of motion in 5D 
theories. We discuss the necessary condition for the large- N limit of the Migdal approach to 
correspond to sensible 5D setups and find which geometries and which boundary conditions 
can be matched. We shall see that the precise manner in which the limit is taken is 
important. 

Finally, we study the relation between Migdal's approach at finite N and deconstruc- 
tion. Deconstruction is a four-dimensional framework that involves a product gauge group 
G N and a set of bifundamental non-linear sigma model fields (the links) [12]. Deconstruc- 
tion models are parametrized by a set of gauge couplings gj and decay constants Vj. For 
large N such setup is related to a 5D gauge theory with the gauge group G, where the 
fifth dimension is latticized. Each choice of gj and Vj on the deconstruction side corre- 
sponds to some 5D warped geometry, and the dictionary between the two frameworks has 
been established [14]. In this paper we quantify the relation between deconstruction and 
the Migdal approach. We show that, given the coefficient of the recursion relation for 
Pade approximants, we can identify the deconstruction parameters that yield the same 
polarization function as the Migdal approach. An unexpected result is that deconstruction 
models directly related to the Migdal approach are non-minimal; they must include the 
kinetic mixing between neighbouring gauge fields from the product group. This, however, 
corresponds to an irrelevant operator in the continuum limit. 

The paper is organized as follows. In Section |21 we present several examples that 
illustrate the connection of the Migdal approach with 5D theories. In Section |3] we review 
the mathematical results connecting Pade approximants of Stieltjes functions to orthogonal 
polynomials. The large-iV limit of Pade approximants of Stieltjes functions is studied in 
Section |3] and the relation to 5D theories is quantified in Section In Section El we discuss 
the connection of Migdal approximation at finite N to 4D deconstruction models. Section 
|7| contains our conclusions, and in the appendix we derive the holographic formula for 
two-point functions in deconstruction. 

2 Migdal's approximation in examples 

For two-point correlation functions, the deep Euclidean limit p 2 — > — oo is a function f(t) 
depending on a variable t = p 2 /fi 2 , where /i is an arbitrary renormalization scale. The 
function f(t) has a branch cut along the positive real axis and contains a perturbative 
piece plus power corrections induced by the condensates. We will be interested in the 
perturbative part only. Migdal proposes to compute Pade approximants of f(t) at some 
Euclidean point — A < 0. The Pade approximant is given by a ratio of two polynomials of 
degrees M and N, such that its Taylor expansion around t = — A matches that of f(t) to 
order (t + X) M+N+1 . Under certain assumptions, this series of approximants converges to 
f(t) when M, N — > oo. The idea of Migdal is to take instead a combined limit N — > oo, 
t — > 0, keeping t = tN 2 and M — N fixed. For finite momentum p, this amounts to sending 
N and \i to infinity with \x — ji/N fixed. We call this the Migdal limit. In this limit, the 
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spacing between the poles of the Pade approximants is controlled by the scale jj, introduced 
in the limiting process. Intuitively, even though the spacing between poles goes to zero at 
large N, so as to reproduce the branch cut, zooming in the small t region simultaneously 
allows to resolve them. It turns out that, for the functions f(t) of interest, the limiting 
expression has only simple poles located at timelike momenta and with negative residues. 
Therefore, the Migdal limit of f(t) can be interpreted as the complete two-point function in 
the large N c limit. By construction, this function has the correct deep Euclidean behaviour. 

In this section we present several examples of the Migdal procedure. The details of the 
calculations are postponed until Section El where we present a generalized approach to this 
kind of computation. 

We start with the example discussed in [10]. The two-point correlation function of two 
conserved vector currents has the general form 

iWp)= (^r^/-) n( P 2 ). (2.i) 

For the vector and axial currents in QCD, the leading perturbative contribution to II(p 2 ) 
at large Euclidean momentum p 2 < is proportional to —p 2 \og(—p 2 /n 2 ). 2 Ignoring mul- 
tiplicative constants, we take f(t + 1) = — log(— t). Next, we approximate f(t + 1) by 
II = t[N/N]f, where [N/N]f = Rn/Sn is the Pade approximant to f(t + 1) at t = — 1 
and Rn, Sn are polynomials in t of degree N. These polynomial can be determined to be 

R N (t) = (t+ 1) N P N , S N (t) = (t+ lfP N , (2.2) 

with P/v the Legendre polynomial of degree N and P/v the associated Legendre polynomial 
(see the next section for the definition of associated orthogonal polynomials). The factors 
(t+l) N in (|2.2j) cancel out in the quotient; their role is only to make Rn and Sn polynomials 
in t. In the Migdal limit, Rn and Sn approach 



R{t) = - log(t) J (2 \ft) + 7tF (2VT) , (2.3) 
S(t) = J (2VI) . (2.4) 



On the other hand, a calculation of the two-point boundary correlator for a gauge field in 
AdSs with Neumann boundary conditions at the IR brane yields [5], in the limit zuv ~~ * 0, 

rrf 2x 2 - log(pzuv) Jo (pzm) + itY (pz m ) 

J (2pzm) 

We see that identifying /1 <-> z^ y , p, <-> 2z^ the result riMigdai = tR/S agrees precisely with 
the corresponding holographic calculation. Note that, according to the dictionary above 
and the definition of p,, 2N corresponds to the inverse warp factor zir/^uv- 

2 We take the principal branch of the logarithm with the branch cut along the real negative axis (ar- 
gument 9 G (— 7r, 7r]), and define accordingly non-integer powers. The signature of the metric is mostly 
minus. Physical amplitudes at time-like momenta are evaluated with p 2 right above the positive real axis. 
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At this point it is natural to wonder why the Pade approximant chooses Neumann con- 
ditions. Our next example shows that we have actually made this choice when identifying 
the function f(i), leaving the factor t outside of the Pade approximation. Indeed, let us 
take instead f(t + 1) = — t log(— t), compute the (N + 1, N) Pade approximant to fit + 1) 
at t — — 1 and define Iljv(t) = [N + l/N]f = Rn/Sn ■ It is clear that this keeps the same 
asymptotic function as before. The reason for increasing the degree of the numerator is 
to improve the convergence at large This technical point will be clarified in the next 
section. The result is 



R N {t) = (t+ 1 



,JV+1 



P 



(1,0) 



S N (t) = {t + l) N P^ 0) 



JV 



1 - t 

T+t 



p. 



(1,0) 



JV 



1-t 



1+t 



(2.6) 
(2.7) 



with (P}? ) P^*'^ (associated) Jacobi polynomials. The Migdal limit yields, up to nor- 
malization, 



R(t) = -ilog(t)-^Ji (2V 7 !) + vrtF!-^ (2Vt 
1 



S{t) = ^J l [2^t). 



(2.8) 
(2.9) 



The result n M i g dai = R/S is identical to the holographic one for large zuv ~" *■ 0, with the 
same identifications as before, but this time with Dirichlet boundary conditions at the IR 
brane. 

Now, let us consider a different asymptotic behaviour in the deep Euclidean regime. 
We assume that the polarization function of two vector currents approaches — (— t) 1 / 2 for 
t — > oo. This behaviour is very different from the one encountered in QCD. Instead, it 
is the prediction of a holographic calculation in an asymptotically flat space. Let us first 
compute Il N (t) = t[N/N) f = tR N /S N with fit + 1) = (-t)~ 1/2 . We obtain 



R N (t) = (t 
S N (t) = (t- 
and in the Migdal limit, 



P. 



-1/2,1/2 



JV 



Afp-1/2,1/2 



lV vp 



JV 



1-t 

T+t 
l - 1 
T+t 



p 



-1/2,1/2 



JV 



1 + t 



(2.10) 
(2.11) 



R(t) = -^sin (2Vt 
yt 



Sit) = cos 2Vt 



(2.12) 
(2.13) 



In this case, IlMigdai = tR/ S agrees with the holographic two-point function in 5D Minkowski 
space with Neumann boundary conditions at the IR brane! Choosing instead f(t+ 1) = 



5 



-t) 1 / 2 and U N (t) = [N/N]f, we arrive at 
R N (t) = (t + l) N 



14 
t 

l + t 



<^ 1/2 (^)+^ /2 - 1/2/ 



(2.14) 



S N (t) = (t + l) N P 1 N /2 >- 1/2 ( 1 —±), (2.15) 



R(t) = -cos (2V t) , (2.16) 
5(f) = -^sin hVt). (2.17) 



with the Migdal limit 



Vi 

As could be already guessed, n Migdal = R/S is in this case the same as the holographic 
two-point function in 5D Minkowski space with Dirichlet boundary conditions on the IR 
brane. 

We have given several examples in which the Migdal procedure in a mysterious way 
reconstructs the full 5D correlation functions from information about the deep Euclidean 
limit. One should be however aware that this "magic" does not always work. For example, 
let us try to reproduce the result for the two-point function of a scalar operator of conformal 
dimension 3, which is dual to a massive scalar in AdS5, with mass M 2 = —3k 2 . The 
asymptotic result of the AdS calculation in the deep Euclidean is, up to a constant term, 
— tlog(— t). This is the same as for gauge bosons, except for the constant term. But the 
Pade approximant is determined essentially by the non-analytic piece of the function, as 
we discuss in the next section. Therefore, choosing f(t + 1) = — log(— t) (and multiplying 
by t at the end) or f(t + 1) = —tlog(-t), we arrive at the same Pade approximants as 
in the first two examples. In the second case, when a local polynomial term is added, we 
reproduce the holographic result for the massive scalar with Dirichlet boundary conditions. 
However, the first choice does not give the holographic function for Neumann boundary 
conditions, but rather one with mixed boundary conditions on the IR brane. 



3 Pade approximants and orthogonal polynomials 

In this section, following mathematical literature [15], we present a more systematical 
approach to Pade approximants. Given a function f(s) with a Taylor expansion at s = 
we can define the Pade approximant as follows. We introduce two polynomials Rfj(s), 
Sn(s) of degrees N + J, N, respectively. We choose them such that their ratio has a 
Taylor expansion at s = that matches the Taylor expansion of f(s) up to terms of order 

S 27V+J+1. 

^f\ = f(s) + 0(s 2N +^) (3.1) 
We also assume S^(0) ^ 0. The Pade approximant is defined as 

[N + J/N] f {s) = ^f\. (3.2) 
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Figure 1: The contour of integration T. 



We will often omit the specification of the function and/or the variable, and write simply 
[N + J/N]. Note that R J N and are determined up to normalization only, and that both 
depend on iV and J. In the following we will restrict the input function / to be a Stieltjes 
function. The input functions we studied in the previous section belong to this class up 
to a certain number of subtractions, as we discuss at the end of this section. A Stieltjes 
function is defined by the Stieltjes integral representation 

m = r (3.3, 

Jo I- su 

where <f)(u) is a bounded, nondecreasing function on < u < oo, with finite real-valued 
moments 



f j= / uPtfyiu), j = 0,1,2,... , (3.4) 
Jo 

and A > 0. An expansion of (J3.3)) in power series at s = gives the Taylor series 

3=0 

which converges in the open disk \s\ < A. We assume in the following that the function 
4>(u) is strictly increasing (and hence the measure strictly positive) except for, at most, a 
discrete number of points. Then, f(s) has a branch cut along A < s < oo. 

There is a remarkable connection between Pade approximants and orthogonal polyno- 
mials, which we derive next (for reviews of orthogonal polynomials see, for instance, [16]). 
From the defining equation (|3.1[) . it is clear that 

d m 

j-^{f{s)S J N {s))\ s=Q = ^ m = N + J + l,...,2N + J. (3.6) 



Then we can use Cauchy's integral formula to write these conditions in the form of contour 
integrals: 

S J N {z)f{z) 



I 



dz m+i = °> m = N + J+l,...,2N + J (3.7) 
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with the path T displayed in Fig. ^ The integral along the small semicircle at A vanishes 
since, from its definition (|3.3|) . f(z) has at most a logarithmic singularity at z = A. The 
integral along the big circle also vanishes when J > — 1. Finally, the remaining integrals 
above and below the branch cut cancel except for the discontinuity in the imaginary part 



of f{z). On the other hand, from ()3.3|) this jump is related to the measure d<f>(u) 
by 



— (f(u + ie) - f(u - ie)) ,\<u<oo. 



Therefore, for J > — 1 and any N > 0, 

-uS J N (u) 



a 



m+2 



0, m = N + J +1. 



2N + J. 



Changing variables to w = u 1 and shifting m — > m — N — J — 1, we arrive at 



dwW J {w)w m {w N S J N {w~ 1 )) =0, m = 0, . . . , JV - 1, 



4>'(u)du 
(3.8) 

(3.9) 
(3.10) 



with W J (w) = w J+1 (p'(w). The factor in parenthesis is a polynomial of degree N. Eq. (J3.1(J|) 
shows that the set of polynomials 



7T 



iV 



(w) = w N S J N (w- 1 ), N = 0,1,... 



(3.11) 



is a system of orthogonal polynomials over the interval (0, A x ) with weight W J (w). We 
shall also use the notation dcj) J (w) = w J+1 d(f>(w) = W J (w)dw. Conversely, 

(3.12) 



S 7T N {S t 



This determines the denominators of the Pade approximant [N + J/N] up to normalization. 
To calculate the numerators, let us define the function 



F J (w) 



and the associated orthogonal polynomials 



d(f) J (u) 



w — u 



w — u 



which have degree N — 1. Then, 



where 



AM 



7T J N (W)F J (W) 



p j nM + KH 



(3.13) 



(3.14) 



(3.15) 



w — u 



1 

w 
w' 



d<p J (u) 

A- 1 



u 
w 



w 



u 
w 



-m—1 



u 



N 



u/w 



1X N 



(3.16) 
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is of order w 1 at large \w\. On the other hand, expanding the integrand in the definition 
of F(s) at u = we find 



J( s~ l ) = s~ J I f(s) - \ ]fjS j ) . (3.17) 




Therefore, we obtain 



m = y: fit 1 + + °( s2N+J ^ ( 3 - 18 ) 



j=0 

so that the numerator of the Pade is 

j 

Rn(s) = S J N (s) £ /,V + ^p^" 1 ). (3.19) 

We see that it is a polynomial of degree N + J, as required. 

Putting all the pieces together, we have shown that the Pade approximants of a Stieltjes 
function can be expressed as 

[N + J/N] f (s) = ^f^ + S -^f^. (3.20) 



3=0 



An important consequence of their relation to orthogonal polynomials is that the Pade 
approximants (with J > — 1) of Stieltjes functions have only simple poles inside the open 
interval (A, oo), and all the residues are negative. This is the basic property which allows 
to make contact with large-iV c and with holography at the classical level. It also follows 
from (J3.16|) that the series of Pade approximants [N + J/N] of f(s) converge to f(s) in 
the limit N — > oo, for all J > —1, in the region \s\ < X. (This is not necessarily so 
for the Pade approximants of a general function.) The rate of convergence is geometric, 
as A(w) ~ O (w - ^" 1 ). This can be extended, with weaker rate of convergence, to the 
domain C\[A, oo). 

The orthogonal polynomials on the real line obey a three-term recurrence relation of 
the form 

ir N+ x(w) = (a N w + b N )n N - c N ir n -i, (3.21) 
with > 0, real constants and cn = (aTv^Ar)/ (azv-i^v-i), where 

h » = M*»)= JdwW{w) ( 3 - 22 ) 

is the squared norm. We shall call ajv, b^, cn = aw/aN-i the coefficients with normal- 
ization Hn = 1. The initial conditions for the recurrence relation are 7r_i = (or cq = 0) 
and 7r = C a nonvanishing constant. From their definition, it is clear that the associated 
polynomials satisfy the same recurrence relation but with initial conditions po = and 
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Pi = (IqC J dwW(w). From (j3.21j) . a three-term recurrence relation for Pade numerators 
and denominators (|3.23|) follows: 

Ti +1 {s) = « + b J N s)T^(s) - c^T^s) , (3.23) 

where T = R, S and we have explicitly indicated the dependence of the coefficients on 
J. The equation is identical for S and R, but the initial conditions are not. These 
follow from the ones for the orthogonal polynomials and associated orthogonal polynomials, 
respectively. For instance, when J = — 1, Rq 1 = 0, R^ 1 = const, Szl = 0, Sq 1 = const. 
Note the factor s 2 in the last term, which makes the equation different from the one for 
orthogonal polynomials. This extra factor is very relevant for our purposes, as it has a 
non-trivial role in the Migdal limit. It can be shown that, because the weight function 
W has a compact support, the combinations l/a J N and b J N /a J N will be bounded. It is 
possible to reverse our line of argument taking us from a Stieltjes function to the recurrence 
relation (j3.23|) for its Pade numerators and denominators. According to Favard's theorem, a 
series of coefficients {a n , b n } with bounded l/a J N and b J N /a J N determines a unique compactly 
supported measure d<fr(u), and a system of orthogonal polynomials with respect to this 
measure, such that these coefficients appear in their recurrence relation. The measure, in 
turn, determines the Stieltjes function. 

All the two-point asymptotic functions f(s) which appear in the conformal approxima- 
tion (which in QCD applies to the leading and subleading perturbative contributions) 
are of the form f n (s) = — (s — A) n log(A — s), with n a positive integer or zero, or 
fu(s) = (— 1)^1 +1 (A — s) u with v a positive non-integer real number and [u] the integer 
closest to v with [u] < v (the entire part for positive v). The variable s is related to 
t = p 2 /fi 2 by s = t + A. These functions are analytic at s — 0, and their Maclaurin 
series have radius of convergence A. On the other hand, the divergent behaviour at large 
\s\ can be taken care of by performing n + 1 or [v] + 1 subtractions, respectively. The 
subtracted functions are Stieltjes functions, and hence their Pade approximants satisfy all 
the properties we have just derived. For instance, for fi(s) = (1 — s) log(l — s), we need 
two subtractions: 

fi(s) = 



So, fi/s 2 is a Stieltjes function, analytic inside the open circle of radius 1 centered at 
s = 0. This equation is nothing but ()3.17)1 with /i(s) = -F J (s _1 ) and J = 2. Therefore, 
the function 

Kn(s) = /i(0) + sf[(0) + s 2 [N - l/N] h (s), (3.25) 

is exactly the same as the Pade approximant [N + l/iVj/^s). It is clear that this is 
generalized to J = m — 1 in the case of m subtractions. 

Consider any f u (s) with real v (possibly integer) and compute the [N + J/N] as above, 
with J > [u] — 1. Changing variables to x = 2w — 1 = (1 — £)/(! + t) we find a weight 



fi(s) - h(0) - sf[(s) 
(1 - s) log(s - 1) + s 




(3.24) 
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W J {x) = (1 - xf{\ - xY J - u *> with x e (-1, 1). Therefore, n J N {x) = > and p J N {x) = 

p(y>J v ) ^ an( j foe Pade approximant is 

[N + J/N]U = + J,j-*),W ■ ( 3 - 26 ) 

i=o ^iv { s ) 

From their integral representation given by (|3.14j) . the associated Jacobi polynomials can be 
written in terms of hypergeometric functions. Note that we should only consider functions 
f v with v > — 1. In fact, when v < — 1, the function /„ is too divergent at s = A to be a 
Stieltjes function. Of course, the Pade approximant can still be calculated, but it will not 
share the good (physical) properties we have derived in this section. 



4 Migdal's limit 

Let us study the Migdal limit of the Pade recurrence relation (|3.23|) . We take A = 1 and 
call r = \/t (with r positive for positive t). the Migdal limit is N — > oo and r — > with 
fixed J and f = tN. To find the limit of equation ()3.23|) . we write (suppressing the index 
J) Tjv(s) = T(N, r), treat # as a continuous variable and expand 

T N±1 (N, r) = T(N, r) ± ^T(iV, r) + ~T(JV, r) + . . . . (4.1) 
Then, keeping only terms up to two derivatives, we get a differential equation of the form 



d 2 _l- CiV (l + r 2 ) 2 9 , 1- fljv-Ml + r 2 ) + c^(l + r 2 ) 2 



9iV 2 1 + c N (l + r 2 ) 2 dN 1 + c n {1+t 



2\2 



T(N,t) = 0. (4.2) 



Let us assume now that the coefficients of the recurrence relation can be expanded at large 
N as 

a N = aW+aV±- + aW-L + ..., (4.3) 

Jy JS Z 

and similarly for and c^. Then, the second order differential equation (|4.2J1 has a finite 
non-trivial Migdal limit if and only if the following conditions are met: 

c<°> = 1, 

2 _ a (°) _ &(o) = o, (4.4) 
C (D _ a (i) _ 6 (i) = o. 

As long the norm of the orthogonal polynomials can also be expanded at large N (with 
a finite number of terms with positive powers of N), the first condition is satisfied. The 
second condition is ensured by Rakhmanov's theorem [18]: if the measure d(f) is supported 
in [—1,1] and > almost everywhere in [—1,1], then it belongs to the Nevai class 3 
with a(°) = 2 and = 0. The assumptions of the theorem are fulfilled by the measure 

3 Measures in the Nevai class are those with finite limits a n — > a,(°>, b n — > [17]. 
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of Stieltjes functions, in the variable x = 2\~ 1 w — 1. Even though the coefficients of 
the recurrence relation in this variable are different, when A = 1 (for which (|4.4j) apply), 
the value of a n + b n is unchanged. On the other hand, 

a (o) = s (o) 6 (o) = 5(0) for the 

class of norms just mentioned. The third condition is more restrictive. Even if it is norm 
dependent, it cannot be adjusted without spoiling our assumption that the coefficients can 
be expanded in 1/N. We check explicitly below that this condition is fulfilled by Jacobi 
polynomials. 

If all three conditions are met, we find for N — > oo the following second-order differential 
equation: 



-^-c«^ + (2-&(°)) + (c< 2 ) 
dr z t clt 



,( 2 ) 



In the case of standardized Jacobi polynomials Pi a '^ with any a and (3 we have, going 



i> ,2, )i 



T(f) 



0. 



(4.5) 



back to the variable w, a {0) = 4, a (1) = -2, 6 (0) = -2, b {l) = 1, 



1 and = -1 



so 



we see explicitly that all the conditions f)4.4j) are directly satisfied. In this case, (14.5)1 reads 



d 2 1 d a 2 ' 
— ^ + -^ + 4- — 



dr 



f dr 



T(f) 



0. 



This is a Bessel equation, with general solution 



T(f) 



C 1 J a (2f) + C 2 y Q (2f) 



(4.6) 



(4.7) 



Actually, in the more general case of Eq. ()4.5|) we can write T(f) = f^+c (1) )nv(f) and 
rescale the variable to a = {^/2 — jffi/2j f. Then, V(a) obeys (j4.6|) with f — > a and 

a 2 = a (2) +6 (2) -c (2) + (l + c (1) ) 2 /4. Note that the common factor f ( 1+c(1) )/ 2 will cancel out 
in the quotient R/S. Actually, this factor comes from the normalization of the orthogonal 
polynomials. On the other hand, the rescaling of f amounts to a rescaling of /L Therefore, 
we see that, as far as the Migdal limit is concerned, and if the limit exists, it is sufficient 
to consider the limit of Jacobi polynomials and work with equation (|4.6|) . 

When the third condition is not fulfilled, the recurrence relation does not have a good 
Migdal limit. In these cases, one could still try to find a continuous differential equation 
by modifying the way in which the limit is taken, and this was actually done (in a different 
language) in Migdal's original paper [8]. In the following we consider only the simplest 
case in which the Migdal limit, as defined here, is finite. We have seen that this reduces 
to studying the Pade approximants of "conformal" functions /„ and /„. 



5 The AdS /Migdal correspondence 

In this section, we describe the general relation between the Migdal approximation and 
extra dimensions, for two-point functions. We give a simple argument showing that the 
Migdal limit unavoidably gives a result which corresponds to a holographic calculation 
in certain 5D geometries. The converse result does not always hold: in some cases the 
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holographic results cannot be obtained from a Migdal approximation to their asymptotic 
Euclidean functions. In order to simplify the notation we consider correlators of scalar 
operators, and comment at the end on the generalization to higher spins. 

We start with the field theory (Migdal) side, and show that analyticity 4 and the con- 
dition of finite the Migdal limit, together with information about the deep Euclidean limit 
and the leading infrared behaviour, completely fix the two-point function n Migdal (t). We 
have seen at the end of the previous section that a good Migdal limit necessarily gives 
numerators R and denominators S which satisfy the differential equation (J4.6|) . Therefore, 
R and S have the form (|4.7|) . with coefficients which are functions of r. 5 

Now, let us impose the asymptotic value of the two-point function. The limit |f| — > oo 
is equivalent to N — > oo with fixed r. Hence, for any non-positive r 2 , the Pade approximant 
converges to fir 2 + 1). Therefore, the Migdal limit of the Pade approximant must have 
the form 

R _ f(r 2 + l)J tM (2f)+A(r)H^(2f) 



(5.1) 

with Hjp the first Hankel function, which goes to zero exponentially for large positive 



S J„(2f) +B(t)H { V (2f) 

3I funct 

imaginary part of the argument. 

Next, we impose that R and S be analytic in p 2 at p 2 = 0. This fixes the functions 
A{t) and B{j). Recall that for integer n, J n (z) is an entire function and Y n (z) equals 
2/ it log z plus an entire function, whereas for any /i, z~ tl J tl {z) is entire. Consider first 
f n (s) = — T 2n log(— r 2 ). Then, the index of the Bessel functions must be an integer, and 

R = -r"- m log(r 2 )J w+m (2f) + nr n - m Y n+m {2f) 
S t n m J n+rn (2f) 

On the other hand, for f u {s) = (— 1)'^ +1 (— t 2 ) u and v a non-integer real, analyticity 
requires fi = v + m with integer m, and 

R = J_,_ m (2f) 
S r v m J u+m (2f) 

Finally, we use the fact that the Pade approximants have no poles or zeros at p 2 = 0, 
and assume that this still holds in the Migdal limit. 6 The leading behaviour at small p 2 is 
R/S ~ p~ 2m in all cases. Hence, we see that m = and this gives the final result for R/S. 
Now, let us define the Migdal two-point function as = t~ l [N + J/N] t it, with integer I. 
As discussed above, we should only consider I > — 1 — v. Since R/S has no poles or zeros 
at p 2 = 0, dal will have a zero of degree — I if I < 0, and a pole of degree I if I > 0. On 
the other hand, the asymptotic behaviour is independent of /. So, we can use the result 

4 The numerators and denominators are analytic in s since they are finite limits of polynomials, and 
hence convergent Taylor series at s = 0. 

5 We are treating r = p/ \i and f = p/jl as independent variables. In the following we also need to use 
the behaviour in the variable p 2 , for fixed /x, jl. 

6 This may be proven using the asymptotic distribution of zeros of the orthogonal polynomials. 
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above with m = I and find 



MigdaI ~ T^T) ' (5 - 4) 

11 Migdal - . /a ' 

J I/+ j(2Vt) 

in the integer and non-integer cases, respectively. Conversely, we see that we can repro- 
duce a low-energy behaviour ~ p~ 21 in the two-point function simply by choosing / in the 
definition of the Migdal two-point function. 

We turn now to the holographic calculations in 4+1 dimensions. In order to keep 4D 
Poincare invariance, the geometry must be a warped direct product of Minkowski times a 
one dimensional space /, which can be chosen flat. In order to have a discrete spectrum, 
/ must be compact. As long as the warp factor is strictly monotonic, one can define 
coordinates in which the metric is manifestly conformally flat: 

ds 2 = i{zf (dx^dx^ - dz 2 ) . (5.6) 

In the coordinate z, I = [zuv>Zir]- AdS geometry corresponds to £(z) = (kz)~ l ; in this 
case, the UV (IR) boundaries at zuv ( z ir) hide the AdS boundary (horizon) at z = 
(z = oo). The holographic prescription to calculate correlation functions of field-theory 
operators at large N (and strong t'Hooft coupling) is given by the AdS/CFT correspon- 
dence [2] : calculate the value of the action for on-shell bulk fields with fixed UV boundary 
values, which act as sources for the dual operators; then, differentiate functionally with 
respect to the sources and put them to zero. For two-point functions, the on-shell action 
reduces to a boundary term. For scalars, 



n(p)= IhnJX 



d z G(z,p) 
[ G(z lP ) 



+ counterterms > , (5.7) 



where G(z,p) is the bulk-to-boundary propagator in a mixed position-momentum repre- 
sentation, fulfilling some specified boundary conditions on the IR boundary and free on 
the UV. X is a p-independent factor which cancels a divergent factor in the non-analytic 
part in the limit Zuv ~~ * 0. The remaining counterterms in (|5.7J1 form a polynomial in p 2 
which cancels the poles at zuv = 0. The propagator G satisfies the equation of motion 
for the dual field. We can write a generic IR boundary condition as G(z\^p) = 0, with 
f(z,p) = Ki(p 2 z 2 )f(z,p) + H2(p 2 z 2 )zd z f(z,p). Such a boundary condition can be obtained 
including mass, kinetic and higher- derivative terms localized on the IR boundary. Indeed, 
higher derivatives of G in the variable z can be written in terms of G and d z G using the 
bulk equation of motion. Let J and Y be two independent solutions of the equation of 
motion. Then, we can write 

G(z m ,p) = J(z m ,p)Y(z,p) - Y(z m ,p)J(z,p) . (5.8) 



14 



Taking the limit, the holographic formula in which the IR conditions are manifest has the 
form 

m = AMAz^p) - A 2 ( P )Y( Zm ,p) + bcal termg ^ (5 g) 

Mip)J{zm,p) ~ A 4 (p)Y(z m ,p) 
Our aim is to relate the differential equation satisfied by the limit of the Pade numerator 
and denominator to a differential equation for K(zm,p) in the variable zjr, where K is any 
linear combination of J and Y . Then, meromorphicity, which follows from the discreteness 
of the spectrum in a compact space, will imply that the holographic two-point function 
will be the same as the Migdal one for some value of /. The function / in the Migdal 
approach corresponds to the limit Zir — > oo of (J5.7|) . with p 2 < 0. In fact, for the functions 
we are considering, which only depend on p 2 /n 2 , we should take as well a low-energy limit 
in which all the scales but |p| (and 1/zir, which we have sent to zero) go to infinity. Then, 
the equivalence of the complete function will hold only in this limit (with finite zir). 
In terms of the length variable L = 2//J, Eq. (J4.6|) reads 



d 2 . 1 d 2 a 



H Yv 

dL 2 LdL L 2 



T(Lp/2) = 0. (5.10) 



Note that dimensional analysis plus the fact that p 2 appears only as p 2 T completely fix 
the form of this equation. On the other hand, the equation of motion of a scalar of mass 
M is 

[d 2 + 3 (d z log£(z)) d z + p 2 - e(z)M 2 ] <t>{z,p) = 0. (5.11) 

In the case of Dirichlet boundary conditions, K\ = 1, k 2 = 0, and K(zm,p) = K(zi^,p) 
obeys the same equation as (p(z,p): 

+ 3 (^ir lo ^m)) d zm +p 2 - ^z lK ) 2 M 2 ] K(z m ,p) = 0. (5.12) 

We see that, because we are using conformal coordinates, the normalization of the term 
with p 2 in (|5.12j) is the same as in (J5.1U|) . If M ^ 0, complete agreement with the Migdal 
equation f)5.1Uj) requires that £(z) = (kz) -1 , with k a constant with dimensions of mass. 
Therefore, the space must be a slice of AdS with curvature k. In this case, 3<9 2lR log£ = 
— 3/zir. To go to the normalization of (j5.1Uj) we only need to write K(zir, p) = z 2 R H(zin, p) . 
Then, for Zm = L and m 2 = M 2 /k 2 = a 2 — 4, the equation for H is exactly the same as 
Eq. (|5.1(J|) . From the AdS/CFT relation between the conformal dimension of the operator 
A and the mass of the dual field, we see that A = a + 2, as it should (remember that a 
is the exponent of the asymptotic two-point function, which is determined by conformal 
invariance). On the other hand, in the case M = 0, we find agreement with Migdal 
equation for any £(z) = (kz)* 1 . Then, we define K(zi&,p) = Zjr 3?? H (zir, p) and identify 
a 2 = (3rj — l) 2 /4. Note that when a^0 there are two different values of t] which give the 
same a. Of course, for rj ^ — 1 the geometry is not asymptotically AdS, and the AdS/CFT 
dictionary linking masses to conformal dimensions must be modified. In the particular 
case 7] = we have flat space, and we see that the same Dirichlet two-point function, 
behaving asymptotically like — (— t) 1 / 2 , is found in flat space with a massless scalar and 
in AdS with m 2 = —15/4. At any rate, adjusting the mass parameter m we can always 
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reproduce Migdal's differential equation in AdS with Dirichlet boundary conditions in the 
IR. Moreover, Dirichlet boundary conditions give a correlator which (for scalars) has not 
a zero nor a pole at p 2 = 0. Indeed, we can rescale the field with a ^-dependent factor 
such that the mass term in the equation of motion cancels. This shows that any zero- 
mode must be flat. But then, the Dirichlet condition forces it to vanish. This is true for 
any UV boundary condition. For UV Dirichlet (Neumann) conditions, this is telling us 
that the (inverse) two-point function does not have a pole at zero momentum. 7 Therefore, 
we conclude that IR Dirichlet boundary conditions correspond to I = in the Migdal 
approach. 

Consider now Neumann boundary conditions, K = d z K. As long as M = 0, and for 
any 77, K satisfies the following second-order differential equation: 



zm • m P 4d 



K(z 1R ,p) = 0. (5.13) 



Writing K = z^ 3r, " 2 H, we reproduce Migdal equation (J5.10|) for a 2 = (37/ + l) 2 /4. 

Thus, aDirichlct ~ "Neumann = _3 ?7. For positive a and 77 < -1/3, a D irichlet - "Neumann = 1- 

Furthermore, in this case there would be a zero mode if Neumann boundary conditions were 
also used in the UV, so the two-point function has a simple zero. Therefore, IR Neumann 
boundary conditions in a massless theory corresponds to choosing / = — 1 in the Migdal 
approach. On the other hand, if M 7^ 0, it is possible to write a second-order differential 
equation for K, but with coefficients which are not analytic in p 2 . So, this equation is not 
of the Migdal form and we cannot reproduce the holographic function within the Migdal 
approach. 

Let us study next mixed boundary conditions with Ki = 1 and t% 2 a constant, and 
assume 77 = —1 to start with. The n, the function K satisfies a differential equation of the 
Migdal form if and only if k 2 = 2±v/4 + m2 , After rescaling, we find a 2 = 5 + m 2 =i=2-\/4 + m 2 , 
respectively. It turns out that this boundary condition, which can be understood as arising 
from a mass term localized on the IR brane, is automatic when the scalar field is the 
super symmetric partner of a fermion or a gauge boson [19]. In this fine-tuned 
mixed UV boundary condition arises as well, such that a zero-mode results. Therefore, 
this corresponds again, up to local terms, to a Migdal function with I — — 1. On the other 
hand, if M = and 77 arbitrary, k 2 = 1/(37/ — 1) and a 2 = [3(77 — l)/2] 2 . Note that 
the remaining solution which would correspond to infinite k 2 , is the one studied in the 
Neumann case. 

It is also possible to reproduce any integer value of I with I + ctDirichiet > — 1 by choosing 
adequate analytic functions K\ and k 2 . This can be proven showing that a good differential 
equation for K is obtained for discrete values of the coefficients (depending on the mass) in 
the expansions of K\ and k 2 , and studying their behaviour at small momentum. Describing 
this in detail would be lengthy, so we simply observe that these properties follow quite 
straightforwardly from the differential-recursion relations of Bessel functions and leave the 
details to the interested reader. 



7 For the relation between holographic correlators and connected correlation functions of the AdS theory, 
including the propagator, see the third reference in [5]. 
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So, to summarize, every Migdal approximation of a two-point function in the (finite) 
Migdal limit can be reproduced by an AdS calculation with a given mass and fine-tuned 
boundary conditions. The Euclidean asymptotic behaviour is determined by the mass, 
whereas the different discrete choices of the parameter I, which controls the leading be- 
haviour at p 2 = 0, correspond to different discrete IR boundary conditions in AdS. The 
converse is not true: not for any mass and IR boundary condition can one find an equiva- 
lent Migdal approximation. One can alternatively reproduce the Migdal calculations using 
a different warped geometry when the mass of the dual field vanishes. In all cases, the 
infrared Migdal parameter jl is proportional to the inverse of the position of the IR bound- 
ary in the conformal coordinates (the ones for which the metric is manifestly conformally 
flat). 

All this discussion can be readily extended to higher integer spins. Ultraviolet conformal 
invariance determines the form of the two-point function in the UV as 



where i,j represent Lorentz indices and is a tensor, polynomical in p^. We choose 
n such that Z^jp 2n be adimensional. For instance, for a vector operator of conformal 
dimension A, 



and n — 1. For a conserved current, A = 3 and we get a transverse function. We can 
directly apply Migdal's approximation to /, and define 



This keeps the tensorial form dictated by conformal invariance. On the other hand, the 
holographic calculation will preserve the conformal tensor structure if the geometry is that 
of a slice of AdS. It will also preserve this form for any metric in the case of completely 
antisymmetric tensors which are dual to p-forms, due to gauge invariance. Therefore, 
the results we have obtained for massless scalars can be generalized to higher p-forms, 
and in particular to gauge fields. For instance, the unphysical example in Section |21 with 
asymptotic behaviour H(p 2 ) ~ p corresponds to a gauge field in flat space. An AdS 
calculation with a vector field with m 2 = —3/4 (corresponding to conformal dimension 
5/2) and adequate boundary terms would reproduce the scalar function U(p 2 ), but not a 
transverse tensor. Finally, in extending the discussion for scalars to tensors, one should also 
take into account that the coefficients in the second term of the equation of motion ()5.11|) 
will be different, which leads to a different normalization and to a different relation between 
m and a (or, equivalently, between m and A). For vector bosons, the coefficient 3 is changed 
to 1 and a D irichict = Vl + m 2 . 

6 Deconstruction and holography 

In the previous section we have studied the Migdal limit of Pade approximants. In this 
one, we relate the approximants with finite N to deconstruction models. In particular, 




(5.14) 




(5.15) 




(5.16) 
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this makes explicit the mechanism by which Migdal correlators approach the holographic 
ones, and how a discrete version of the holographic formula (J5.7|) is realized by the Pade 
approximants. For definiteness, we will stick to the case of gauge bosons — dual to conserved 
currents — and discuss Dirichlet and Neumann conditions only. 

We consider a deconstruction model corresponding to the moose diagram sketched 
in fig. 13 It involves a chain of N SU {Np) groups with the gauge fields = Ajj a T a , 
j = 1 . . . N. The groups communicate with nearest neighbours via bifundamental non- 
linear sigma model field Uj, referred to as the links. They are unitary Np x Np matrices 
with determinant equal to 1. The role of the last non-linear sigma field U is to control the 
boundary conditions for the gauge field (the analogue of IR boundary conditions in 5D). 

At the left end of the chain we singled out the "boundary" gauge field A®. Similarly as 
in AdS / CFT, this boundary field is interpreted as an external current probing the dynamics 
of the "bulk" model. Here, the bulk refers to the remaining gauge fields A^ , j > 1. The 
latter will show up as resonances in the boundary field correlation function. 

In the following we calculate the two-point correlation function of the boundary fields in 
deconstruction. We do it first in the more familiar minimal deconstruction set-up, and then 
in what we call tilted deconstruction, which contains additional interactions between the 
neighbouring sites. Correlators obtained in tilted deconstruction turn out to be directly 
related to Migdal's Pade approximants and we work out a dictionary between the two 
approaches. 



6.1 Minimal deconstruction 

The gauge transformations u>j, j = . . . N, act as Aj^ — > UjA^uo^ — id^UjUj' and Uj — > 
ujj^iUjUJj^ . The simplest non-trivial action that is invariant under these transformations 
can be written as 

S = j d^fj(-^tr{i^ + J d'xtriv'D.UD^}, (6.1) 

with DfjUj = dpUj — iAi~ l Uj + iUjA^ , D^U = d^U — iA^U. This action is minimal in 
the sense that the only interactions between various gauge fields come from the covariant 
derivatives acting on the links. The relation of this deconstruction setup to 5D gauge 
theories can be worked out analogously as in [14] . Consider the 5D action for a gauge field 




Figure 2: The moose diagram for our deconstruction setup 
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in the background ds 2 = a 2 (z)dx 2 — b 2 (z)dz 2 : 

(6.2) 

Latticizing the 5th coordinate, z — > z 3 - = jA, d$f(z) — > (f(zj) — f(Zj-i))/ A we obtain: 

S 5D - /^E (-^ trF i(^) + J^) tr ^(^ - M^-i)) 2 ) ■ (6-3) 

This can be mapped onto the deconstruction action (jfj.ljl (in the unitary gauge Uj = 1). 
The warp factors and the lattice spacing translate into the parameters of the deconstruction 
lagrangian according to the following dictionary: 

a (Zj)->— — > K z j)^— , A ^ , 9 5 ^—- (6-4) 

vigj 9j givi vi 

We have fixed a{z{) = b(z\) = 1. In passing we note that discretization in Poincare 
coordinates, b(z) = 1, corresponds to gj = g, while discretization in conformal coordinates, 
b(z) = a(z), corresponds to gjVj = gv. 

We can integrate out all the bulk gauge fields and obtain an effective action for the 
boundary field. At tree-level, the integrating-out amounts to 1) solving the equations of 
motion for the bulk fields in the presence of a background boundary field and 2) inserting 
the solution into the deconstruction action. The details of this procedure are given in the 
appendix. At the quadratic level, the effective action in the momentum space has the form 

Seff = J ^-A A » (-^ + V ~f) KipMp 2 ) . (6.5) 
The polarization function is given by a compact expression 



n(p 2 ) 9 2 v 2 



p 2 - 1 . (6.6) 



where F 3 N is a solution to the equation of motion 

v] +1 + v 2 - 4 W - v*F^ - v 2 +1 F 3 N +1 = , (6.7) 
9 j / 

subject to a boundary condition at j = N. The boundary condition is controlled by the 
parameter v in the lagrangian. The limit v — > leads to a deconstructed analogue of the 
Neumann boundary condition, 

F% +1 = F$ , (6.8) 
while v — > oo corresponds to the Dirichlet boundary condition, 

7*N 



F£=0. (6.9) 



19 



Intermediate values of v correspond to mixed boundary conditions. On the other hand, 
for g —> oo we recover the case of non-dynamical boundary fields, which can be regarded 
as sources. 

There are many apparent similarities between polarization functions obtained in min- 
imal deconstruction and those derived using the Migdal procedure. Let us point them 
out. 



The correlation function is represented as a ratio of two polynomials in p 2 

Rn( P 2 ) 
Sn(p 2 ) 

where 



n(? 2 ) = l4§. (e-io) 



Rn — + ( —2~2 ~ 1 ) $n — Fn (6-H) 
\9o v i J 

Indeed, from the equation of motion ()6.7|) F 3 N is a polynomial of degree N — j in p 2 , 
once we set F$ = const (in the Dirichlet case F N ~ X = const and F 3 N has the degree 
N-j-l). 

2. The numerator and the denominator satisfy a second-order recurrence relation in the 
degree N, and the equation is the same for both. For example, in the Dirichlet case 
it is given by (T = R, S): 



L + ) T N (p 2 ) - ^T^p 2 ) . (6.12) 

v n+i 9n v n+i/ v n+i 



This follows from the fact that the solution F 3 N satisfying the Dirichlet boundary 
conditions can be written as F 3 N = Y^Jj — J^Yj, where Jj and Yj are any two 
independent solutions to eq. (|fi.7j) . 

3. In the limit — > oo and for \p 2 \ v\ and p 2 < 0, the polarization function obtained 
in deconstruction approximates the non-analytic behaviour of the polarization func- 
tion in the corresponding 5D model in the deep Euclidean regime. For example, in 
the deconstructed AdS models one obtains [13] H{p 2 ) ~ — p 2 log(— p 2 /vf), while in 
the deconstructed flat models we find H(p 2 ) ~ p 2 (—p 2 /vf)~ 1 / 2 . 

In spite of these similarities it is not possible to find a precise mapping between the 
Migdal approximation and minimal deconstruction. The reason is that the recurrence rela- 
tions ()6.12j) and (|3.23j) are incompatible. Indeed, the form of the recurrence relation (J6.12JI 
implies that T^{p 2 ) is an orthogonal polynomial in the variable p 2 . On the other hand, 
the numerators and denominators obtained by the Migdal procedure, although related to 
orthogonal polynomials by eq. ()3.12|) . are themselves not orthogonal polynomials. 

In the following we explore a modified deconstruction framework that allows for a 
mapping of the boundary correlators to those obtained using the Migdal approximation, 
lators to those obtained using the Migdal approximation. 
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6.2 Tilted deconstruction 



We modify the minimal deconstruction action (|6.1|) adding a kinetic mixing between 
neighbouring gauge fields, 

S= fd*x Ef =1 (-^tr{F^J + tr^f/^t/j}) + / d±xtr{v*D„U N+1 D^ N+1 } 

~ I d*x Ef =1 ftr{/^'r ; Z-^r; + h.c.} + / d*x-^tr{F»F»} (6.13) 

Such deconstruction setup is also related to a latticized 5D gauge theory in the warped 
background. The difference with the minimal case is that the 5D action must contain a 
higher derivative term breaking the 5D Lorentz invariance: 

S 5 d - / d*x f dz (-^tri£ + ^-Md,A,f - f^trF,„d 2 F, u ) . (6.14) 
J Jo V ^95 95°\ z ) z 9h J 

The dictionary between deconstruction and 5D is now given by: 



A^^S, g 2 -> . (6.15) 

We see that tilted deconstruction corresponds to a 5D setup with some specific higher- 
derivative terms. Nevertheless, the extra term is irrelevant at low energies and a standard 
lowest-order 5D action is recovered. Hence, the effect of the mixing term amounts to a 
renormalization of the coefficients in this action. Viewed as a 4D field-theoretical model, 
tilted deconstruction is healthy as long as the mixing coefficients aj are not too large (so 
that there are no ghosts). 

As in the minimal setup, it makes sense to integrate out the resonances A* j > 1 and 



calculate the effective action for A ^. The polarization function is given by (see | Appendix A 
for a derivation) 



W) V J 9 2 o< 

The last two terms have a form of a local polynomial in p 2 . By adding higher derivative 
terms for the boundary fields to the deconstruction action we could, in fact, obtain an 
arbitrary polynomial in p 2 . 

In tilted deconstruction, F 3 N solve a modified equation of motion 

+ "1 " J) n - («J + FIT 1 - («J„ + ) FP 1 = , (6.17) 

subject to the boundary condition 

(fs + i ~v 2 + j 2 V 2 ) F» = (v 2 N+1 + ^±y) . (6.18) 

In the limit v — > oo we obtain Dirichlet boundary conditions, F$ = 0, while setting v — 0, 
OTV+i/gAr+i = l/^ 2 leads to Neumann boundary conditions, F^ +1 = F^. 

We will prove that, for certain choices of the coefficients aj, the polarization functions 
obtained in this setup are directly related to those obtained by Migdal approximation. 
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6.3 Migdal— deconstruction map 

Let us make the following ansatz for the mixing coefficients: 

-T = ^> 6 - 19 

where \i is an arbitrary scale. We also introduce a new variable, s — 1 + p 2 /fi 2 . The 
equation of motion now becomes 

+ + ^ - 4) - v]sF^(s) - v 2 sF£\s) = . (6.20) 



We find it convenient to discuss the Dirichlet and the Neumann case separately. 

Dirichlet boundary conditions 

We investigate the solutions to eq. (J6.20J) subject to the boundary condition F$ = 0. 
If we set -F/v" 1 = const then s- 7 F^~^~ 1 (s) is a polynomial in s of degree j. We define the 
polynomials 

RM = s^F^s), SUs) = s N F° N+1 (s). (6.21) 

It follows that Sn has degree N, while Rn has degree N+l. From eq. (|6.16J) the polarization 
function can be written as 

n(«) = |44 + /o + A* ( 6 - 22 ) 

and has the form of a Pade approximant with J — 1. Indeed, the numerator and the 
denominator as defined in eq. (j6.21|) satisfy the second order difference equation 

V u N+2 ilN+l u N+2 Un+I u N+2 J u N+2 

subject to the boundary conditions 

Rl(s) = 0, R\(s) = constant, Sl^s) = 0, Sfo) = constant (6.24) 



Equations ()6.23|) and ()6.24|1 follow simply from the fact that the solution F 3 N+1 satisfying 
Dirichlet boundary conditions can be written as F 3 N+1 = Vat+i J? — J^ + iYj, where Jj and Yj 
are any two independent solutions to eq. (J6.20|) . They have exactly the same form as the 
recurrence equations and the boundary conditions for the denominator and the numerator 
in Pade approximation. 

In fact, the recurrence relation in Pade approximation contains three sets of coefficients 
On, &tv and cjy, which define the related orthogonal polynomial. On the deconstruction 
side we dispose only of two sets: and v^. However, since rescaling of the numerator 
and the denominator by the same function does not change the polarization function, 
Migdal approximation and deconstruction are equivalent if the recurrence equations can 
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be brought to the same form after rescaling by an arbitrary function, — > hj^T^. 
Taking this into account leads to the norm-independent consistency conditions: 

gAT+lC^T+2 + v N+l) ^ a N 

/i 2 b N 
9n9n+i v n+i c n 



(6 ' 25) 

Thus, given a^, b]y and C]y defining the orthogonal polynomial corresponding to the Pade 
approximants, we are able to reconstruct the (tilted) deconstruction model that would give 
exactly the same polarization function. Furthermore, using the dictionary (j6.15|) we can 
find the continuum background. 

Let us now investigate what deconstruction model corresponds to the Pade approxi- 
mants found in section 2. Those examples where all associated with Jacobi polynomials 
P^{2s — 1), whose recurrence relation involves the coefficients 

(2N + a + (3 + l)(2iV + a + (3 + 2) 



a N 



>N 



(N + l)(N + a + f3+l) 
(2N + a + (3 + l)(a 2 - (3 2 - (2N + a + (3){2N + a + (3 + 2)) 
2{N + 1)(N + a + + 1){2N + a + (3) 



( i y + a )( iV + / 3)(2iV + a + /3 + 2) 
° N (AT + l)(iV + « + /3 + l)(2iV + « + /5) • { - ) 



For large N the consistency equations can be approximated by: 
9n+M+2 + v n+i) i , a 2 -(3 2 1 



/i 2 2 N 2 

9Wn+iV% + i 1 1-4/3 2 1 
/i 4 4 16 N 2 

They can be approximately solved by 

v 2 N+1 l±2a 



0(1/N 3 ) , 
0{1/N 3 ). (6.27) 



N 



N 



0(l/N 



2\ 



9n+i 1 _ l±2a 
9n N 



0(1/N 



2\ 



9j f - H 1 -^ ^)- (6 - 28) 

The corresponding deconstruction background, at lowest order, does not depend on (3. We 
see that a = 1/2 can be matched to flat deconstruction with g N and independent of 
N. On the other hand, a = 1 is equivalent to deconstruction with 

t&««i^, 9l*\glN, (6.29) 

which by eq. ()6.15|) is deconstruction of AdS 5 gauge theories latticized in conformal coor- 
dinates. 
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Neumann boundary conditions 

The results for the Neumann boundary conditions F^j +l = F$ can be obtained in an 
analogous way, and below we simply give our results. We define the polynomials: 

R%(s) = -^ l s N (sF 1 N (s) - F° N (s)) S%(s) = s N F° N (s) (6.30) 

Both i?jv and Sn are polynomials of degree N in s (one can show that the factor l/(s — 1) 
always cancels out). The polarization function can be written as 

n <«> = 5 (ft + f °) (631) 

and has the form of a Pade approximant with J = 0. The polynomials satisfy a recurrence 
relation that is different from the Dirichlet case: 

rp ( \ f V N+l(9N + 9N+l) + A* 2 A* 2 \ rp / 2\ 9n v N 2rp i \ tao^\ 
T N+1 (S) = 2 2 2 2 S T "(P ) 2 S T N-l(s) ■ (6-32) 

V 9n+i v n+2 9n+i v n+2 J 9n+\V n+2 

In obtaining this equation, the correlation between the coefficients of (|6.2U|) is crucial. 
Adding a mass term would spoil this correlation and the resulting recurrence relation in 
the Neumann case would not be of the Pade form. This agrees with our discussion of the 
continuum extra dimensions. The consistency conditions are given in this case by: 



v2 N+l (g 2 N+i + 9%) 

a 2 b N 



/i 2 b 

^a^I+i = cn_ (6 33) 

/i 4 b N b N ^ V ; 

For Jacobi polynomials on the Migdal side the large- N approximate solution is given by 

_ 1 l±2a , 0{1/N2) 



'N 



N 



% = 1 + + 0(1/N 2 ) 

9% N 

9>n _ l( 1 + ii^ + 0(1/iV 2 ) ] 



/i 2 2 V 2N 

We can see that a = —1/2 is reproduced by flat deconstruction with g^ and independent 
of N, while a = corresponds again to deconstruction of AdS 5 in conformal coordinates, 

v 2 N ^vf^, g 2 N ^ l -glN. (6.35) 

These results are very welcomed, as they show that a single deconstruction setting is able 
to reproduce the values / = and I = 1 in the Migdal approximation, just by a change of 
the IR boundary condition. 
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7 Conclusions 



In this paper we have studied the relations between three different methods of computing 
correlation functions in strongly-coupled large- N c theories: the Migdal approach via Pade 
approximants, the 5D holographic approach via boundary correlators, and the deconstruc- 
tion approach via external field correlators. We have made explicit the connection between 
the Migdal approximation and the other two methods. The key feature of the Pade approx- 
imant is that its denominator and numerator can be expressed in terms of some orthogonal 
polynomial and their associated orthogonal polynomial, respectively. This ensures physical 
properties of the Pade approximants analogous to those of large- N c theories. Furthermore, 
the recurrence relations satisfied by the orthogonal polynomials provide a link with local 
equations of motion in the other two approaches. 

The equivalence between Migdal and deconstruction correlators for finite N gives a 
nice explicit realization of the so-called UV/IR relation between the energy scale and the 
radial position in holographic models. Indeed, in Migdal's approach the successive Pade 
approximants allow an extrapolation of the UV result to lower and lower energies. This 
corresponds in deconstruction to the addition of new sites and links, which in turn generate 
an extra dimension. Another interesting common feature is that a discrete spectrum is 
obtained thanks to a violation of quark-hadron duality, which is introduced by hand, 
either by keeping a constant t = N 2 p 2 /fi 2 in the Migdal limit or by cutting off the space 
with the IR brane. Locality in N implies that this violation is ~ l/t N , and exponentially 
suppressed in the continuum. 

We have considered the large- N Migdal limit, in which N 2 p 2 /fi 2 is constant. This 
forced us to restrict the input functions to the ones that have a conformal form. On 
the other hand, conformality in the UV is related to asymptotically-AdSs spaces. The 
non-trivial fact that we have explored here is that Migdal's approximation extrapolates 
this conformal/AdS character all the way down from the UV to the IR, up to an abrupt 
IR cutoff/brane. 8 This is related to the particular limit we are considering. It would 
be interesting to investigate non-conformal input functions, which arise at higher orders of 
perturbation theory. This would require a different Migdal limit, and it could be speculated 
that a softer IR cutoff would be generated (possibly involving an infinite extra dimension, 
as in [20]). Conformality is also broken by power corrections which, in this context, were 
discussed in [10]. It would also be interesting to study the correspondence for higher-point 
correlators [21]. 

The relation between Pade approximation and holographic calculations in 5D or decon- 
struction could be regarded as a mere mathematical curiosity. However, we expect it to 
have physical consequences as well. In fact, the Migdal program relies on dispersion rela- 
tions and is similar in spirit to the SVZ sum rules, which have a solid theoretical basis. The 
connection with 5D models might shed some light on the unexpected success of AdS/QCD 
models [6]. Furthermore, Pade approximation is often employed as a unitarization method 
to extrapolate the predictions of chiral perturbation theory in QCD [22] and in no-Higgs 

8 We have also seen that a similar extrapolation is at work in other conformally flat spaces in the 
massless case. 
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models of electroweak breaking [23] to higher energies. This approach is complementary to 
Migdal's, for it goes from low to high energies rather than the other way round. 9 On the 
other hand, alternative unitarization procedures using the notion of extra dimensions have 
been introduced more recently: the so-called higgsless electroweak breaking [26] and its 
deconstructed version [27]. Our results suggest that these seemingly unrelated approaches 
could be equivalent, although to prove it we should study the Pade approximants with a 
low-energy, rather than high-energy, input. 
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Appendix A Derivation of the boundary effective ac- 
tion in deconstruction 

We derive here the holographic formula for the two-point correlation function of the 
"boundary" fields in deconstruction. We work in the tilted deconstruction framework; 
the minimal deconstruction result can be obtained by setting ctj — 0. 
The tilted deconstruction action can be rewritten as 




(A.l) 



where the kinetic operator is defined as 

\9j 9j 9j+i J 

+V»v ({v] + v 2 j+1 )5j, k - vftj^k - v* +1 5 j+1>k ) . (A.2) 

For D^, vq = is understood. Our objective is to obtain an effective action for the 
boundary field A ^ after integrating out all resonances with j > 1. At tree-level, this 
can be achieved by solving the equations of motion for the resonances with the boundary 
background field switched on, 

N 

E^" = ° 3>l, (A.3) 

k=0 

9 In [24], it has been pointed out that, because the exact (large- N c ) two-point function is a Stieltjes 
function, if the chiral contributions to all orders were known, then the large- N limit of the corresponding 
Pade approximants would exactly reproduce the full two-point function, at arbitrary momenta. In par- 
ticular, this puts restrictions on the allowed chiral coefficients, directly related to the ones which can be 
derived from dispersion relations [25]. 
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and inserting the solution back into the action (|A.lj) . The solution can be written as 

where F° N (p 2 ) solves the recurrence relation in j, 

+ v) - g Ff, - («J + ^) - (v] +1 + ?^±l) *T = , (A.5) 

subject to appropriate boundary conditions: F$ = in the Dirichlet case and F^ +l = F$ 
in the Neumann case. Inserting the solution back into the action we obtain 

If d 4 p 2 (f i PuPA .4Qn^ 2 V4° i PnP» aO A o F n( q ) ~ F n(Q) \ (a ^ 

Se *-2j J^ v i\{-^ + ^) A » n (P) A » + ^r A » A » — fW) J ■ ( } 

The second term vanishes in the Neumann case, while in the Dirichlet case it is cancelled 
by tree-level exchange of a massless physical Goldstone boson. Finally, the polarization 
operator is given by 

n(/) = fjg (i + (a.7) 



9tn j gm 
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